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Abstract Local regularity of axially symmetric solutions to the Navier- 
Stokes equations is studied. It is shown that under certain natural assump- 
tions there are no singularities of Type I. 
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1 Introduction 

In this paper we will consider local regularity properties of axi-symmetric 
solutions of the 3D Navier-Stokes equations 

d t v + v ■ Vf + Vg — Av = , , 

diw = o . t 1 - 1 ; 

Most of the known regularity theory for these equations (and, in fact, for 
many other equations) is based on optimal estimates for the linear part and 
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on treating the nonlinearity as a perturbation which is (locally) small in a 
suitable sense. An important role in formulating suitable smallness condi- 
tions is played by certain (local) scale-invariant quantities. These are the 
quantities which are invariant under the scaling symmetry v(x,t), q(x,t) — > 
Xv(Xx, \ 2 t), \ 2 q(\x, \ 2 t). The reason why the regularity criteria should be 
formulated in terms of the scale-invariant quantities is simple: The class of 
regular solutions is invariant under the scaling and therefore sufficient condi- 
tions for membership in this class should ideally be also invariant under the 
scaling, or at least they should scale in the correct way, in that the quan- 
tity controlling regularity should not decrease if we scale the solution with 
A>> 1. 

To write down examples of the scale-invariant quantities we recall the 
following standard notation. The points of the space-time R n x R will be 
denoted by z = (x,t). For x £ R n we denote by B(xq,R) the ball {x : 
\x — xo\ < R} and for z = (xo,to) G I" x M we denote by Q(zo,R) the 
parabolic ball B(xo,R)x]t — R 2 ,t [. Here are some examples of the scale- 
invariant quantities for n = 3: 

(1.2 



(1.3 
(1.4 

(1.5 

(1.6 

(1.7 
(1.8 

A typical local regularity result says that, under some natural technical 
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assumptions [J, a point z is a regular point of the solution v if a suitable scale 
invariant quantity X(z , R; v) of the type in the examples above is sufficiently 
small for all R e]0, Rq[. In fact, X can be any of the quantities above with 
the exception of (11.51) . in which case the validity of the corresponding result 



is open. 

At the time of this writing, there is no known scale-invariant quan- 
tity for which an a-priori estimate would be known for general 3D solu- 
tions. In fact, all the known estimates can be traced back to the energy 
estimate, which gives bounds in quantities such as J \v(x, t)\ 2 dx or 

B(x ,R) 

J \Vv\ 2 dz, which do not have the scaling needed for the existing lo- 

Q{z ,R) 

cal regularity theory. This is often quoted as the main stumbling block in 
our understanding of the Navier-Stokes regularity. This statement is prob- 
ably correct, at least as a first approximation. However, even if we assume 
that scale-invariant estimates of natural quantities are available, in many 
cases we are still unable to prove regularity, unless an additional smallness 
condition is imposed. For the quantities of the type X(z , R; v) listed above 
one can show that X(zq,R;v) < C for some C > (not necessarily small) 
implies regularity for (jl.2p and (11.31) . but in the remaining cases the known 
theory requires an additional smallness condition (and, as remarked above, 
the situation with (11.51) is even worse). Moreover, even the proofs of the cases 
(11.21) and (II. 3p rest on the fact that the assumptions imply that a certain 
quantity becomes smallO 



1 Such assumption must include some control of pressure, such as q G L,3. Fortunately, 
such control is available from energy estimates in most situations. 

2 The reason for the difficulties with (II. 5| is that the space-time L2 norm of v is not 
sufficiently strong to control the energy flux (unless one can come up with some surprising 
new property of the equations). Roughly speaking, the energy flux is controlled by the 
I/3 norm of v. Since the energy estimate gives the control of the Lw norm of u, there is 
some gain in regularity and it is natural to try to bootstrap it and try control the energy 
flux by some L r norm with r < 3. This does work, but the borderline exponent r for this 
argument is r = 5/2, still quite far from r = 2 which would be needed for a local regularity 
result with R~ 3 J \v\ 2 dz. See for example [25] . 
QOo.-R) 

3 The finiteness of (|1.2|l implies lim J \v\ 5 dz = 0, which gives us a small 

quantity. The finiteness of Q1.3P implies (for the solutions of the equation) 
lim J \v(x,to)\ 3 dx — 0. This again gives a small quantity, but in this case it is 

R ^°B(x ,R) 

not easy to exploit it, since we essentially have to show that some regularity propagates 



In this paper we study local regularity results for axi-symmetric solu- 
tions of the 3D Navier-Stokes under an assumption that a weakened version 
of quantity (11.71) or, respectively, (jl.8p is finite (but not necessarily small). 
These studies can be thought of as a continuation of the work started in [3] , 
[TO] , and [I]. The exact assumption which we will use to replace fll.8p . in the 
axi-symmetric situation, with the a^-axis as the axis of symmetry, is 



ess sup y xf + x\ \v{x,t)\ < +00 (1.9) 

for some R > 0, where zq lies on the X3-axis and we denote by v(x, t) the 
projection of the velocity vector v (x, t) into the plane passing through x and 
the axis of symmetry 23. Similarly, the exact assumption which will replace 
(11.71) in the axi-symmetric situation, with the X3-axis as the axis of symmetry, 

is 

ess sup y/to — t \v(x, t) I < +00 (1-10) 

(x,t)eQ(z ,R) 

for some R > 0, where z$ and v are as above. Our main results are as follows. 

Theorem 1.1. Assume that v G L^{Q{zq,R)) is an axially symmetric weak 
solution to the Navier-Stokes equations in Q(zo, R) such that there exists an 
associated pressure field q G Lz (Q(z , R)). If, in addition, v satisfies /J1.10\) . 
then zq is a regular point of v. 

Theorem 1.2. Assume that v G L^(Q(zo,R)) is an axially symmetric weak 
solution to the Navier-Stokes equations in Q(zq,R) such that there exists 
an associated pressure field q G La(Q(zQ,R)). Suppose that v is essentially 
bounded in the space-time cylinders ofthefromB(xo,R)x]to — R 2 ,t'[for each 
t' < t , where the bound may depend on t' . If, in addition, v satisfies U.9\) . 
then Zq is a regular point of v. 

These are local versions of the main results in the paper [ID] . Similar 
(but not identical) results also appeared in [3] and 0]. 

For completeness, we formulate another theorem, which is a local version 
of the corresponding global regularity result in [IT] and [30J . 

Theorem 1.3. Assume that v G L^(Q(zo,R)) is an axially symmetric weak 
solution to the Navier-Stokes equations in Q(zo,R) such that there exists 



backwards in time. 



an associated pressure field q G L3(Q(zq,R)). Suppose that v is essentially 
bounded in the space-time cylinders of the from B(xq, R)x]t — R 2 ,t'[ for each 
t' < to, where the bound may depend on t' . If, in addition, the field v has no 
swirl, i. e. v = v, then zq is a regular point of v. 

On a conceptual level our method will be close to the one used in [10], and 
will rely on the Liouville-type theorems established in that paper. However, 
certain important technical parts will be treated in a different way. 

We first recall some terminology related to the Liouville-type results for 
the Navier-Stokes proved in [10]. An ancient solution of the Navier-Stokes 
equation is a solution defined in R"x] - oo, 0[. We are interested in ancient 
solutions with bounded velocity, see Definition [231 Non-zero solutions of this 
form can be generated by a natural re-scaling and limiting procedures at a 
potential singularity, see Section 2. The definition of the ancient solutions 
still allows for the "parasitic solutions" of the form u(x,t) = b(t) (for any 
bounded b: ] — oo,0[^ M. n ), with the corresponding pressure p given by 
p(x, t) = —b'(t)-x, see Remark 12.61 To exclude these solutions (which - under 
some natural assumptions - cannot arise from the re-scaling procedures, see 
Theorem 12.81) we introduce the notion of the ancient mild solutions. These 
are the ancient solutions which satisfy the natural representation formula 

u{t) = S(t - t )u(to) + / (S(t - s) Pdiv (u(s) ® u(sj) ds (1.11) 

Jt 

for some sequence of times to — * — oo, where S is the solution operator for 
the heat equation and P is the Helmholtz projection onto the div-free fields. 
(We remark that the usual integration by parts shows that the integral on 
the left-hand side is well defined for u G Loo-) See [TQ] for details. The 
strongest conjecture regarding the Liouville-type results one can make about 
the Navier-Stokes equations is the following: 

Conjecture (L): The velocity Geld of any bounded mild ancient solution of 
the Navier-Stokes equations is constant. 

The conjecture was proved for n = 2 and also for axi-symmetric solutions 
in 3D, provided the additional decay condition 



'xl + xl\v{x,t)\ <C inM 3 x]-oo,0[ (1.12) 

is satisfied, see [10J. 



What would be the implications of the validity of Conjecture (L) for 
the regularity theory? Roughly speaking, if Conjecture (L) is valid, then all 
the problems discussed above concerning regularity in the presence of a scale- 
invariant estimates are solved. Indeed, the re-scaling procedure preserves any 
scale-invariant estimate, and typically the estimate will also be preserved in 
the limiting process. Therefore as a result of the re-scaling we get, in the 
limit, a non-zero bounded mild ancient solution for which a scale invariant 
quantity is finite. Conjecture (L) would leave only one candidate for the mild 
ancient solution - namely a non-zero constant velocity field. However, this 
possibility is typically not compatible with a finite scale-invariant bound. 

We can summarize the above as follows: 

scale invariant „, . , T . . ., 

. + Conjecture (L) =>- regularity. 

Singularities for which some scale-invariant quantity is bounded are often 
called Type I singularities. (The most common definition of Type I singu- 
larities uses quantity (11.71) .) 

While we do not really know what the likelihood of Conjecture (L) being 
true is for the general 3D solutions, we are quite confident that the conjecture 
is indeed true for the axi-symmetric solutions. The axi-symmetric case of 
Conjecture (L) would imply much stronger results than Theorems II. II and ll. 21 
above. However, we have not been able to fully prove Conjecture (L) in the 
axi-symmetric case so far. 

Our method of proof of Theorems 11.11 and 11.21 is can be described as 
follows. Roughly speaking, we will show that, on the solutions of the equa- 
tions, the assumed scale invariant bounds imply that all the other important 
scale-invariant quantities are bounded, and these bounds, together with the 
known partial regularity theory ([Tj [T2], [T7]), lead relatively easily to the 
bounds required by the Liouville theorems in [10]. 

The idea that a bound of one scale-invariant quantity should lead (for the 
solutions of the equations) to bounds on other scale invariant quantities is of 
course not new. However, examples from some other elliptic/parabolic PDEs 
show that these issues can be subtle. For example, in the theory of harmonic 
mappings or the harmonic map heat flow we do have a scale- invariant a-priori 
bound, which corresponds to a bound of quantity (jl.5p . However, it is known 
that singularities can still arise, and therefore the bound corresponding to 
(11.21) (which is known to imply regularity in that situation) cannot be derived 
from (the analogue of) (II. 5p . 

6 



We now informally explain the main steps of the proof. One is that the 
swirl component of the velocity field, v v = v ■ e^ satisfies a scalar parabolic 
equation which enables one to gain some regularity. To explain this, we need 
to introduce the following simple notation. Let e\, e 2 , e 3 be an orthogonal 
basis of the Cartesian coordinates x\, x 2 , £3 an d e g , e v , e 3 be an orthogonal 
basis of the cylindrical coordinates g, <p, x 3 chosen so that 

e e = cos (pei + sin (pe 2 , e v = — sin ipei + cos (pe 2 , e 3 = e 3 . 

Then, for any vector-valued field v, we have representations 

v = v^ = v 1 e 1 + v 2 e 2 + v 3 e 3 = v Q e Q + v v e v + v 3 e 3 . 

Next, letting / = gv v , we have 

d t f + v-Vf = Af-- 1 f. (1.13) 

Q OQ 

We would like to prove a Loo-bound on /. Such a bound will give us enough 
information about v — v so that, oversimplifying slightly, we can replace v by 
v in our assumptions. The L^ bound for ( 11.131) does not follow from general 
parabolic theory, since the general theory requires more regularity than we 
have. However, it is known that if the drift term in equations such as ( 11.131) 
is div-free, one can prove the L^ estimate for / with weaker assumptions 
on the coefficients. See for example [6] for the elliptic case and [32] for the 
parabolic case. 

Another important step in the proof is conceptually the same as deriving 
an estimate for the quantities (J1.4J1 and (11.61) from the boundedness of ( II. 7J) . 
This can be done by bootstrapping the energy inequality. This idea was used 
for example in [23] . The technical details are somewhat complicated, but the 
main idea can be explain at a heuristic level as follows. To simplify notation, 
we will use Q(R) for Q(0,R) = Q((0,0),R), Q for Q(l), B(R) for B(0,R), 
and B for B(l). 

We first note that (11.71) implies a bound on the (scale-invariant) quantity 
i? 1_ T _ ~||f || s ,z,Q(fl) for / < 2 and s > 1. Here, || • || Si /,q(_r) is the norm of the 
mixed Lebesgue space L s j(Q(R)) = Li(—R 2 ,0;L s (B(R))). 

Let [w]]n(m) denote the parabolic energy norm in Q(R), i. e. 

( u Iq(R) = ess su Pie]-R 2 ,0[ll' u ('' i )llL2(B(R)) + I|V m IIl 2 (Q(R)) • 



To avoid technicalities, let us pretend that the pressure satisfies \q\ ~ 
|t>| 2 . In reality it is not quite true and in the rigorous proof one has to deal 
with this, but the procedure is well understood. Therefore our simplifying 
assumption |g| ~ \v | 2 is reasonable for the heuristics. We will now work with 
R — 1, but we can scale the calculations to any R > 0, if we divide all the 
involved quantities by the powers of R which make them scale- invariant. 

The local energy inequality implies 



v 



\q ~ IMIl 3 (Q(2)) + IML 3 (Q(2)) • I 1 - 14 ) 



We can now "bootstrap" this inequality. There are some technical com- 
plications coming from the fact that we have Q(2) on the right-hand side 
of (11.141) but only Q on the left-hand side. Such problems come up of- 
ten in local regularity theory of elliptic and parabolic equations, and it is 
quite well-understood how to deal with them by suitable iteration proce- 
dures. Therefore, cheating slightly, we can pretend that we actually have Q 
on the right-hand side of inequality fjl.14)) : 

(1.15) 
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(1.16) 

for suitable a, (3 > 0, a + (3 — 1, and use this in inequality (11.151) . We see 
that when a < 2/3, we can estimate [v\q in terms of the norm \\v\\ sl q. (The 
process also works for a = 2/3 provided \\v\\ sl is sufficiently small.) 

It remains to determine the correct exponents a, (3 in (11.161) . Denoting 
by 2* the Sobolev exponent of the space W 1,2 (i. e. 1/2* = 1/2 — l/n), we 
have by the Holder inequality 

Hull < Il7;ll ai ll7;H a2 Hull" 3 (1 17^1 

1 1 ^ 1 1 3,3 — H f; ll2,ooll f; ll2*,2ll f; lls,Z > \ L - Lt ) 

where ai,«2,«3 are non-negative numbers satisfying 

«i + a 2 + a 3 = 1 

Ql , Q2 , 03 = 1 (1 18) 

2 ~ 2* ~ s 3 l- L -- LO /' 

02 I «3. _ I 

2 ' I 3 



By Sobolev imbedding, we have from (11.171) 



^llwo)<MS 1+Q2 H<tQ> (1-19) 



and we see that (jl.lfip holds true with a = a,\ + «2, (3 = otz- Therefore the set 
of the parameters s, I for which the iteration procedure works is given by the 
condition that equations fll. 18j) for 0^1,02,03 have a non-negative solution 
with as > 1/3. Solving (11.181) for n = 3, we obtain 
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One can check easily that the conditions a% > 0, «2 > 0, 0:3 > 1/3 are 
equivalent to 



1+1 > 2 

s ^ I — 3 ' 

f + T > 1, (1-20) 

! + ! < 2 - 

In the plane with coordinates x = - and y — \, the last set of equations 
describes a thin triangle contained in the first quadrant. It is easy to see 
that one can choose a suitable I < 2 and s > 1 for which these conditions are 
satisfied. 

Above we worked with parabolic balls of radius of order 1. It is clear 
that the calculations can be "scaled" to Q(R) and Q(2R) if we divide the 
L s ;-norms by suitable powers of R to obtain scale invariant quantities. It is 
also clear that of the restrictions in fll.20p . only the last one is crucial, since 
by using Holder inequality one can always move to higher exponents, as long 
as the power in the correct scaling factor remains positive. 

This finishes our explanation of the heuristics behind the second step 
of the proof. We did not mention one more complication. Since our main 
assumption involves only v and the information about v — v is obtained from 
equation (11. 13)) . we have to use different one set of parameters /, s for the v 
component of v and a different set for v — v. However, this is a technical 
issue which does not change the heuristics. 

Once we know that the scaled energy-type quantities are bounded, it is 
not difficult to derive the bounds which we need in the version of the Liouville 
conjecture for axi-symmetric solutions which was proved in [10j . 

9 



Another aim of the paper is to give an alternative approach to certain 
technical issues arising in the study of bounded ancient solutions. In the ap- 
proach here, we do not use the exact representation formulae which were used 
in jlO] . It turned out to be quite convenient to describe differentiability prop- 
erties of bounded ancient solutions in terms of certain "uniform" Lebesgue 
and Sobolev spaces, compare with [T5]. We hope that both approaches are 
of interest, and complement each other. 

2 Preliminaries 

In this Section, we recall known definitions of (weak) solutions to the Navier- 
Stokes. 

Definition 2.1. A weak solution to the Navier-Stokes equations in a domain 
O C M n x]ti,t 2 [ is a divergence free vector-valued field v G L^ioAG) satisfying 



(v ■ d t w + v®v:Vw + v- Aw)dx dt = 
o 

for any solenoidal vector-valued field w e Cq°(0). 

An important family of weak solutions is given by v(x,t) = Vh(x,t) 
where h : O — > K satisfies Ah = 0. (Dependence on t can be arbitrary). 
This example shows that further assumptions are needed to obtain some 
regularity of solutions in the time direction. 

Very often, we shall study local regularity of solutions to the Navier- 
Stokes equations in the unit parabolic ball Q = Bx] — 1,0[, where B = 
B(l) = B(0, 1). It is not a loss of generality because of the Navier-Stokes 
scaling. 

In local analysis, the most reasonable object to study is so-called suitable 
weak solutions, introduced by Caffarelli-Kohn-Nirenberg in pQ . We are going 
to use a slightly simpler definition of F.-H. Lin in [T7] 

Definition 2.2. The pair v and q is called a suitable weak solutions to the 
Navier-Stokes equations in Q if the following conditions are satisfied: 

veL 2tOC (Q)nW^(Q), geL f (Q); 

v andq satisfy the Navier-Stokes equtions in the sense of distributions; 

10 



for a. a. t G] — 1,0[ ; the local energy inequality 

t t 

/ ip(x,t)\v(x,t)\ 2 dx + 2 I I ip\Vv\ 2 dxdt' < I I i\v | 2 (A^ + d t <p)+ 



-1 C -1 c 



+v ■ Vtp(\v\ 2 + 2q)\dxdt' 

holds for all non-negative cut-off functions if G C^°(]R 3 x M.) vanishing in a 
neighborhood of the parabolic boundary ofQ. 

Here, the following functional spaces have been used: 
L a<l (Q) = Li(-1, 0;L s {B)l L s = L 8>8 , 
Wlf{Q) = {v, Vv G L S ,(Q)}, W^° = WA°, 

Kt(Q) = i v > Vu ' ^ V 2 v G L Sil (Q)}, Wf 1 = WH 

The norm of the space L s> i(Q) is denoted by || • \\ s ,i,q- 

For further discussions of Definition 12. 2\ we refer the reader to papers [T2] 
and [23j. 

In what follows, we shall assume v and q satisfy the following standing 
conditions: 

pair v G L 3 (Q) and q G Ls(Q) satisfies the Navier-Stokes equations 

in the sense of distributions; (2.1) 

veL oc {Bx]-l,-a 2 [) Voe]o,i[; (2.2) 

there is a number < r\ < 1 such that v G L^Qx), (2.3) 

where Qi = Bix] — 1,0[, B x = {r\ < \x\ < 1}. To explain why there is no 
loss of generality, we first notice that the pair v and q, satisfying conditions 
f)2.1l) - fl2.3l) . is in fact a suitable weak solution to the Navier-Stokes equations 
in Q. It is certainly true in Bx] — 1, — a 2 [ but condition ( 12.11) allows us to 
extend this property to the whole cylinder Q. 

Consider now an arbitrary suitable weak solution v and q in Q. Let 
S C -Bx] — 1, 0] be a set of singular points of v. It is closed in Q. As it 
was shown in pQ, V 1 (S) = 0, where V 1 is the one-dimensional parabolic 
Hausdorff measure. Let us assume that 5^0. We can choose number Ri 
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and R 2 satisfying < R t < R 2 < 1 such that S n Q{R{) \ Q(R 2 ) = and 
S n B{R 2 ) x] - R%, 0] ^ 0. We put 

£ = inf{£ : (x,t) E S(l B(R 2 )x] - R\, 0]}. 

Clearly, (x ,t ) G 5 for some x G B(R 2 ). In a sense, t is the first sin- 
gular time of our suitable weak solution v and q in Q(Ri). Next, the one- 
dimensional Hausdorff measure of the set 

S to = {x* G B(R 2 ) : (£*, to) is a singular point} 

is zero as well. Therefore, given xq G St , we can find sufficiently small 
< r < A/-R2 + *o such that 5(i ,r) <s S(i2 2 ) and dB(x ,r) n S^ = 0. 
Since the velocity field t> is Holder continuous at regular points, we can ensure 
that all conditions of type (I2.1l) -( j2~3l hold in the parabolic ball Q(zq, r) with 
z o = (^o, *o)- We may shift and re-scale our solution if xo = and r^ 1. 

In our investigations of regularity of suitable weak solutions, the partic- 
ular case of weak solutions, see Definition 12.11 plays a crucial role. Here is 
the corresponding definition. 

Definition 2.3. (Jffll ) A bounded divergence free field u G L OQ (Q_;W a ) is 
called a weak bounded ancient solution (or simply bounded ancient solution) 
to the Navier-Stokes equations if 

/ (u ■ d t w + u ®u : Vw + u ■ Aw)dz = 
Q- 

o 

for any w G Co°(Q-)- 

Here, we have used the following notation: 

Q_ = R n x]-oo,0[, CZ > (Q-)={veC?(Q..), divw = 0}. 

The notion of bounded ancient solutions is not quite satisfactory for our 
purposes since it allows "parasitic solutions" u(x, t) = b(t), where b is an 
Loo-function. (These correspond to harmonic function h(x,t) = b(t) ■ x.) 

The important subclass of bounded ancient solution was introduced in 
[TO] . It consists of the so-called mild bounded ancient solutions, i.e., bounded 
ancient solution satisfying the representation formula (II. lip . In [10J, we also 
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showed that one has a natural decomposition u(x,t) = w(x,t) + b(t), where 
w is given by the right hand side of (II. lip on a suitable interval ]£i,£*[- (In 
particular, w is Holder continuous.) 

In this paper, we give another proof of the decomposition of arbitrary 
bounded ancient solutions into regular and singular parts, see Section 5. It is 
based on recovery of a pressure field associated with a given bounded ancient 
solution. To formulate our theorem about the pressure, we need to introduce 
certain functional spaces 

By L m (fl) and W^fi), we denote the usual Lebesgue and Sobolev spaces 
of functions defined on Q e M n . We also need parabolic versions of Sobolev's 
spaces: 

W^°(Q T ) = {\v\ + \Vv\eL m (Q T )}, 

W 2 m \Q T ) = {\v\ + \Vv\ + \d t v\ + \V 2 v\ E L m (Q T )} 

where Qt = f2x]0,T[. The norm of the space L m (fl) is denoted by || ■ || m ,n- 
We also going to use the following "uniform" spaces (compare with |15j): 

Cm(Q-) = {||/|U m (Q_) = sup ||/|| m ,Q( 20 ,i) < +00}, 
V\&°(Q_) = {||/llwfc°(Q_) = SU P ll/ll^°(o(«,i)) < + °°>' 

z a eQ- 

v^; 1 (g-) = {||/IUy(Q_)= su p WfWwZiHQM) < +0 °>- 

zo£<5- 

To define the regular part of the pressure, we recall the known fact, see 
e.g. [28]. Given F = Loo(IR n ; M nx?1 ), there exists a unique function pp G 
BMOiW 1 ) such that [pf]b(i) = ([/]n is the mean value of a function / over 
a spatial domain Q e ~R n ) and 

Ap F = -div div F = F ijAj in M 3 

in the sense of distributions. Moreover, function pp meets the estimate 

||Pf||bmo(m™) < c(n)||F|| 00)K n. 

So, given a bounded ancient solution u, we define a regular part of the pres- 
sure Pu®u- 
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Theorem 2.4. Let u be an arbitrary bounded ancient solution. For any 
number m > 1, 

|Vu| + \V 2 u\ + \Vp u(g>u \ G C m (Q-)- 

In addition, for each t < 0, there exists a function b to G I><x>(£o ~~ l^o) with 
the following property 

SUP IKIUoo(to-Mo) < c ( n ) < +°°- 
<o<0 

If we let u to (x,t) = u(x,t) + b to (t), (x,t) G Q to = R n x]t - l,t [, then, 
for any number m > 1, 

SU P \\ ut[, \\w? n \Q(z ,l)) < C K n ) < +°°- 

z o =(xo,to),x o m n ,t o <0 

Moreover, for each t < 0, functions u and u to obey the system of equations 

dtu to + div u® u — A-u = — Vp u( g, u , div u = 

a. e in Q to . 

Remark 2.5. The first equation of the latter system can be reduced to the 
form 

d t u + div u®u — Au = —Vp u ^ u — b' to in Q to , 

which is understood in the sense of distributions, b' to (t) = db to (t)/dt. So, the 
real pressure field in Q to is the distribution p u ® u + b' tQ ■ x. 

Remark 2.6. We can find a measurable vector-valued function b defined on 
] — oo,0[ and having the following property. For any t < 0, there exists a 
constant vector ct such that 

sup||6-Ct ||L oo (t -i,t ) < +°°- 

t <0 

Moreover, the Navier-Stokes system takes the form 

d t u + div-u (g) u — Au = — V(p u( g, u + b' ■ x), divu = 
in Q- in the sense of distributions. 
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Remark 2.7. Bounded ancient solutions with b' = were introduced in J72J/ 
and called mild ancient solutions. They were systematically studied in [W^ 
and in particular it was shown there that mild ancient solutions are infinitely 
smooth. 

Our interest in bounded ancient solutions comes from their appearance 
as natural limits of suitable re-scaling procedures at potential singularities. 
In the context of solutions to the Navier-Stokes equations in R n x]ti,t2[, this 
was studied in [10J. 

We shall now show that re-scaling procedure also works for potential 
singularities of local suitable weak solutions. The important point will be 
that, even in the local (in space) situation, the solutions arising from the 
re-scaling procedure at a potential singularity are still mild bounded ancient 
solutions. To be more precise, let us consider local solutions of the Navier- 
Stokes equations satisfying assumptions (I2.1I) - (I2.3I) and introduce functions 

G{t)= max \v(x,t)\, M(t) = sup G(t). 

cceB(ri) -l<r<t 

Assume that there are singular points of v which are located somewhere on 
the set {(x, 0) : \x\ < r{\. By known regularity criteria (see e.g. [25]), we 
have 

M{t) > 



-t 

for some e > and thus 

M(t) -> +oo 

if t — > 0— . We can construct a sequence tk such that tk 6] — 1, 0[, tk < t k +i, 
tk — > 0, and 

M(t k ) = G(t k ) = \v(x k ,t k )\ -> +oo 

for some x& G B{r\). 

Next, we scale v and q the following way 

u k (y, s) = X k v(x, t), p k (y, s) = X 2 k q(x, t), 

where x = x k + X k y, t = t k + A^s, and A^ = 1/M k . The ball B{r) is mapped 
by the change of variables onto B(—x h /X k , r/X k ) and if r £]r 2 , 1] then, given 
any R > 0, 

k 

B(J*)cB(-f,f 

\ Afc Aa 
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for sufficiently large k. For scaled functions u k and p k , we know 
u k and p k satisfy the Navier-Stokes equations and 



» A |<! inZ?(-^-,^]x 



l-t fc 
A! 



.0 



(2.4) 



|u fc (0,0)| = l. (2.5) 

Theorem 2.8. For eac/i a > 0, the sequence {u k } is uniformly Holder con- 
tinuous on the closure ofQ(a) for sufficiently large k and a subsequence {u kj } 
of {u k } converges uniformly on compact subsets ofR n x] — oo, 0] to a mild 
bounded ancient solution u with |w(0,0)| = 1. 

Remark 2.9. The main point of the theorem is that the limit u is a mild so- 
lution, i.e., the parasitic solutions cannot appear in this re-scaling procedure. 

Proof of Theorem 12 .8I Our solution v and q has good properties inside 
Qi. Let us enumerate them. Let Q2 = B 2 x] — r|,0[, where < Ti < 1, 
B2 = {n < T2 < |x| < 02 < 1}. Then, for any natural k, 

z = (x,t) 1— > V k v{z) is Holder continuous in Q 2 ; 

gGL f (-r 2 2 ,0;C fe (g 2 )). 

The corresponding norms are estimated by constants depending on ||f ||3,q, 
ll<?IU,Q5 IM|oo,Qi) an d numbers k, r 1; r 2 , a 2 , r 2 . In particular, we have 


max / \Vq(x,t)\idt < c x < 00. (2.6) 

x€B 2 J 



Proof of this statements can be done by induction and founded in [5], [12], 
and [IB]. 

Now, let us decompose the pressure q = qi + g 2 . For q x , we have 



Aqi(x, t) = — divdiv xb(x)v(x, t) ®v(x,i) 



xew, -1 < t < 0, 



where xb(x) = 1 if x & B and Xb(x) = if x ^ B. Obviously, the estimate 


/ / \qi(x,t)\2dxdt < c / \v\ 3 dz 



Q 
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holds which is a starting point for local regularity of q\. Using essentially the 
same bootstrap arguments, we can show 

o 

max / \V qx(x,t)\ 2 dt <c 2 <oo, (2.7) 

x£B 3 J 

-I 
where B = {r 2 < r 3 < \x\ < a 3 < a 2 }. From ( 12.61) and (12. 7ft . it follows that 

o 
max / \Vq 2 {x,t)\^dt <c 3 <oo. (2.8) 

x£B z J 



But clearly q 2 is a harmonic function in B, thus, by the maximum principle, 

we have 

o 

max / \Vq 2 (x, t)\?dt < c 3 < oo, (2.9) 

es(r 4 ) J 



xeB(r 4 ) 



where r 4 = (r 3 + a 3 )/2. 

Let us re-scale each part of the pressure separately, i.e., 

rf!(v,s) = \lqi(x,t), z=l,2, 

so that p fc = p k + p 2 . As it follows from ( 12.9f) . for p|, we have 

o 
sup / |V^(2/, S )|irf S < c 3 a|. (2.10) 

y£B(-x k /\ k ,r 4 /\ k ) J 

-(rl-t k )/X 2 k 

The first component of the pressure satisfies the equation 

AyPiiVi s ) = ~ div y div y (XB(-x"/x k ,i/x k )(y)u k (y, s) <8> w fc (i/, s)), J/ G M 3 , 
for all possible values of s. For such a function, we have the standard estimate 

lbi(-,s)|| BMO (K3)<c (2.11) 



for all a 6] - (1 - t k )/Xl, 0[. 

We slightly change pf and £>2 setting 

Pi(y, s) = Pi(y, s) - [Pi]s(i)(s) ^(y, s) = p k 2 (y, s) - \p 2 ] B (i)(s) 
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so that \pi\b(i)(s) = and [p 2 ]b(i)(s) = 0. 

Now, we pick up an arbitrary positive number a and fix it. Then from 
(I2.10p and (j2.1ip it follows that for sufficiently large k we have 



/ \p k \^de + I \plfide < c 4 (c 2 ,c 3 ,a). 

Q(a) Q{a) 

Using the same bootstrap arguments, we can show that the following estimate 
is valid: 

\\ uk \\c-{Q{a/2) < C 5 (c2,C 3 ,C4,a) 

for some positive number a < 1/3. Indeed, the norm ||w fc |lca(Q( a /2)) is es- 
timated with the help of norms ||w fc ||L 00 (Q(a))) and ||p fc ||L 3 (Q(a)), where p k = 

Pi + p\- Hence, using the diagonal Cantor procedure, we can select subse- 
quences such that for some positive a and for any positive a 

u k -^u in C a {Q{a)), 
Pi^Pi, inIrs(Q(o)), \Pi]b(i)(s)=0, 

p\ ^P 2 in L l(Q(a)), [P2\b(i)(s) = 0. 

Moreover, u is a bounded ancient solution with the total pressure p = Pi+p 2 , 

where p x = p u<S)U . 

Next, for sufficiently large k, we get from (j2.10p that 



sup / \Vp k 2 {y,s)\Us<c z \l. 

PR(a) J 



Q(a) 

Hence, Vp 2 = in Q(a) for any a > 0. So, p 2 (y, s) is identically zero. This 
allows us to conclude that the pair u and p u ® u is a solution to the Navier- 
Stokes equations in the sense of distributions and thus u is a nontrivial mild 
bounded ancient solution satisfying the condition |w(0,0)| = 1. Theorem 12.81 
is proved. 

3 Axially Symmetric Suitable Weak 
Solutions 

Without loss of generality, the problem of local regularity of weak solutions 
(not necessary being axially symmetric) to the Navier-Stokes equations can 
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be formulated as follows. Let us consider a pair of functions v G L 3 (Q) 
and q G La(Q), defined in the unit space-time cylinder Q = Cx] — 1,0[, 
where C = {x = (x',x 3 ), x' = (xi,x 2 ), \x'\ < 1, \x 3 \ < 1} is the unit spatial 
cylinder of M 3 , which satisfies the Navier-Stokes system in Q in the sense 
of distributions. The question we are interested in is under what additional 
conditions on v and q, the space-time origin z — (x, t) — is a regular point 
of v. By the definition, the velocity v is regular at the point z — if there 
exists a positive number r such that v is essentially bounded in the space-time 
cylinder Q(r). Here, Q(r) = C(r)x] — r 2 ,0[ and C(r) = {\x'\ < r, \x 3 \ < r}. 
In contrast to traditional setting, we replace the usual balls with cylinders, 
which is quite convenient in the case of axial symmetry. As usual, we set 

v = v g e e + v 3 e 3 v = v^e^ 

for v = v e e Q + v v e v + v 3 e 3 . 

We reformulate our main results for these canonical domains. The general 
case is obtained by re-scaling. 

Theorem 3.1. Assume that functions v G L 3 (Q) and q G Ls(Q) are an 
axially symmetric weak solution to the Navier-Stokes equations in Q. Let, in 
addition, for some positive constant C , 

\v(x,t)\<-j= (3.1) 

for almost all points z = (x, t) G Q. Then z = is a regular point of v. 

Theorem 3.2. Assume that functions v G L 3 (Q) and q G La(Q) are an 

axially symmetric weak solution to the Navier-Stokes equations in Q. Let, in 
addition, 

ueLoo(Cx]-l,-a 2 [) (3.2) 

for each < a < 1 and 

\v(x,t)\<^- (3.3) 

\x'\ 

for almost all points z — (x,t) G Q with some positive constant C. Then 
z = is a regular point of v . 

It is well-known due to Caffarelli-Kohn-Nirenberg that if z = (x, t) is 
singular (i.e., not regular) point of v, then there must be x' = 0. In other 
words, all singular points must belong to the axis of symmetry which is axis 
x 3 . 
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Lemma 3.3. Assume that functions v G L 3 (Q) and q G Lz(Q) are an 
axially symmetric weak solution to the Navier-Stokes equations in Q. Let, in 
addition, condition A3.2\) hold. Then following estimate is valid: 

/_3_ 
Igv^l^dz) , (3.4) 

0(3/4) 

where 

M=( [ |U|*dz)* + l. 

0(3/4) 

For the reader convenience, we put the proof of Lemma 13.31 in Appendix 
II, see also [2] and [3]. Here, we would like to notice the following. 

Remark 3.4. Under the assumptions of Lemma \3.3\ the pair v and q forms 
a suitable weak solution to the Navier-Stokes equations in Q. Hence, the right 
hand side of (3.4 ) is bounded from above. 

We recall that the Navier-Stokes equations are invariant with respect to 
the following scaling: 

u(x, t) = Xv (Ax, A 2 t), p(x, t) = \ 2 q(\x, X 2 t) 

So, new functions u and p satisfy the Navier-Stokes equations in a suitable 
domain. 

With some additional notation 

C(x ,R) = {x G 1R 3 || x = (x'jXz), x' = (xi,x 2 ), 

\x'-x' \ < R, \x 3 -x 03 \ < R}, C(R)=C(0,R), C = C(1); 
z=(x,t), z = (x ,t ), Q(zo,R) =C{x ,R)x]t - R 2 ,t [, 
Q(R) = Q(Q,R), Q = Q(1), 
we introduce certain scale-invariant functionals: 

A(zq, r; v ) = ess sup - / \v (x, t)\ 2 dx, 

t -r 2 <t<t r J 

C(xo,r) 

E(z ,r;v) = - / \Vv\ 2 dz, D(z ,r;q) = — / \q\*dz, 

Q(zo,r) Q(zo,r) 
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C(z ,r;v) = — / \v\ 3 dz, H(z ,r;v) = — / \v\ 2 dz, 

Q(zo,r) Q(*o,r) 

M Si i{z ,r;v) = — / ( / M s dx) 3 d£, 

t -r 2 C(xo,r) 

where k = l(- + | — 1) and s > 1, 1 > 1. As it was shown in |25j . the following 
inequality holds 

C(z , r; /) < cA"(zo, r; /)(M.,,(zb, r; f))™(E(z , r; /)+ 

+ ff(zb,r;/)) 2 £ 1 , (3.5) 

where 

" = mt + T- 2 )' m = a U + T-2)' 

provided 

3 2 3 (I 1 1 l-i , . 

- + --->max ---,--- . 3.6 

s / 2 1 2 s s 6 J 

Actually, inequality (13.51) is but the result of application of Holder's inequality 
and special Galiardo-Nireberg's inequality. 

The essential technical part of the proof of Theorem 13.11 is the following 

lemma. 

Lemma 3.5. Under assumptions of Theorem \3.1\ we have the estimate 

A(z b , r; v) + E(z b , r; v) + C(z b , r; v) + D(z b , r; g) < d < +oo (3.7) 

for all Zb and for all r satisfying conditions 

z b = (6e 3 , 0), beR, \b\ < -, < r < -. (3.8) 

.4 constant C\ depends only on the constant C in A3. 1\) . |M|l 3 (Q), o,nd 
IMI-MQ)- 

Proof By Lemma [3.31 and by Remark 13.41 we have two inequalities: 

A(0, 3/4; v) + E(0, 3/4; v) < C 2 < +oo, (3.9) 

W\\v<p(x,t)\ < C 2 for &.&.z= (x,t) E Q(l/2). (3.10) 
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Constant C<i depends on the same arguments as constant C\. 

It follows from A3.5J) . that, for s\ = % and l\ = 10, inequality (13. 10J) takes 

the form 

58 1 

mi = — , /ii = — , M s Zl (2; 6 ,r;?) < cC™, 

7 58 
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C(z b ,r;w) < cA™(z b ,r;v)(C™)ss(E(z b ,r;v) + A(z h ,r;v))» (3.11) 

provided conditions ( 13.81) hold. 

To treat v which is the other part of the velocity v, we chose numbers 
S2 = 4 and li = y. Then, for the same reasons as above, we find 

!0 3 ,, , , « 

™2 = y , M2 = j^, M S2ii2 (z b , r;v) <cG 7 , 

C{z b , r- v) < cA^(z b , r- v)(C^)^(E(z b , r; v) + A(z b , r; v)) iff (3.12) 

for all z b and r satisfying conditions ( 13.81) . 
Adding (IXTTD and (I3TT2|1 . we show 



C(z&, r; u) < c C(z 6 , r; v) + C(z 6 , r; u) I < 



<c^A-8(z b ,r;v)C> 9 (E(z b ,r;v) + A(z b ,r;v))-s+ (3.13) 

+AT3 (z 6 , r; v)C§ (S(^ 6 , r; v) + A(z b , r; u))wj 

for the same 2& and r as above. 

Applying Young's inequality in (13.131) . we arrive at the important estimate 

C(z b , r; v) < e(E(z b , r; u) + A(z b , r; «)) + f^e, C, C 2 ), (3.14) 

provided conditions (I3.8P hold. In (13.141) . the positive number e is a param- 
eter to pick up later. The rest of the proof is routine. In addition to (I3.14J) . 
we consider the local energy inequality 



E(z b , r/2; v)+A(z b , r/2; u) < c(C* (z b , r; w)+C(z fe , r; u)+D(^, r; g)J (3.15) 
and the decay estimate for the pressure field 

D(z b ,g;q)<c -D(z b ,r;q) + (-) C(z b ,r;v) . (3.16) 

Lr \g/ J 
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Here, z b and r satisfy conditions (13.81) and < g < r. If we let 
£(r) = E(z b , r; v) + A(z b , r; v) + D(z b , r; q), 



then, for a fixed small positive number $, one can derive from (I3.15P and 
(13. 161) the following estimate 

£(tfr) < cfd (z b , 2$r; v) + C(z 6 , 2$r; v) + £>(z 6 , 2tfr; g) + 



< c 



+tf J D(z b , r; g) + ^C{z b , r; v)J < 

1 1 2 ■ 

$£>(z 6 , r; g) + -^C(z b , r; v) + -rCa (z 6 , r; v) 



tf2 > - ' ^ 

Now, the last two terms on the right hand side of the latter inequality can 
be majorized with the help of (j3.14p . As a result, we have 

E{<dr) < c(# + ^)s(r) + f 2 (e,#,C,C 2 ). 

We first chose $ so that c& < j, pick up e to provide the inequality J§ < j, 
and then we find 

£{tir)<\s{r)+h{C,C 2 ). 

The latter inequality can be easily iterated. After simple calculations, we 
derive the relation 

E(z b , r; v) + A(z b , r; v) + D{z b , r; q) < c(a(0, 1/2; v) + E(0, 1/2; v)+ 

+D(0,l/2;q) + f 3 (C,C 2 ; 

with z b and r satisfying conditions (I3.8p . Lemma [3.51 is proved. 

To prove Theorem 13. 2\ we need an analogue of Lemma 13.51 Here, it is. 



Lemma 3.6. Under assumptions of Theorem \3.2\. estimate ( 3.1 ) is valid as 
well with constant C\ depending only on the constant C in 113.3]) , \\v\\l 3 (Q), 

and ll<?IU 3 (Q)- 
2 
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Lemma 13.61 is proved in the same way as Lemma 13.51 and even easier 
because main inequality (I3.14p can be established with the help of the case 
s — S\, I — l\ only. 

As it follows from conditions of Theorem 13.21 and the statement of Lemma 
13.31 the module of the velocity field grows not faster than C/\x'\ as \x'\ — > 0. 
Moreover, the corresponding estimate is uniform in time. However, it turns 
out to be true under conditions of Theorem 13.11 as well. More precisely, we 
have the following. 

Proposition 3.7. Assume that all conditions of Theorem \3. 1\ hold. Then 

\v(x,t)\<^ (3.17) 

\x | 

for all z = (x,t) e Q(l/8). A constant C\ depends only on the constant C 
in (E2]j ; |M|l 3 (Q), and \\q\\ Lz {Q)- 



Proof In view of (13. 5p . we can argue essentially as in 
Let us fix a point xq G C(l/8) and put r = \x' \, bo = £o3- So, we have 
r$ < I and |6q| < |- Further, we introduce the following cylinders: 

'Pro = ( r o < W\ < 2r , |x 3 | < r }, V^ = {r /4 < \x'\ < 3r , |x 3 | < 2r }. 

K ( & o) = K + & oe 3 , K (b ) = V 2 ro + b e 3 , 
Ql (b )=V} (b )x]-rlO[, Q 2 rQ (b )=V? (b )x] - (2r ) 2 ,0[. 
Now, let us scale our functions so that 

x = r Q y + b e 3 , t = r%s, u(y,s) = r v(x,t), p{y,s) = r%q(x,t). 

As it was shown in [26], there exists a continuous nondecreasing function 

$ : M + -f R + , R + = {s > 0}, such that 

sup \u(u, s)\ + |Vw(y, s)\ < $( sup / \u(y, s)\ 2 dy 

P?(0) 

+ / \Vu\ 2 dyds+ I \u\ 3 dyds+ / \p\?dyds\ (3.18) 

Q?(o) Q?(o) Q?(o) 
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After making inverse scaling in (I3.18p . we find 

sup r \u(x,t)\ + rl\\7u(x,t)\ < $( cA(z ho , 3r ; v) + cE(z bo , 3r ; v)+ 



^Ql {b ) 



+cC(z bo ,3r ;v) + cD(z bo ,3r ;q)^j < $^4cC] 

It remains to apply Lemma 13.51 and complete the proof of the proposition. 
Proposition 13.71 is proved. 



4 Proof of Theorems 13.11 and 13.2 



Using Lemmata 13. 3[ 13.51 13.61 Remark 13. 4[ Proposition 13.71 and scaling argu- 
ments, we may assume (without loss of generality) that our solution v and q 
have the following properties: 

sup (A(0, r; v) + E(0, r; v) + C(0, r; v) + D(0, r; qj) = A 1 < +oo, (4.1) 

0<r<l V / 

ess sup \x'\\v (x, t)\ = A 2 < +oo. (4.2) 

z=(x,t)eQ 

We may also assume that the function v is Holder continuous in the comple- 
tion of the set Cx] — 1, — a 2 [ for any < a < 1. 
Introducing functions 

H(t) — sup \v(x, t)\, h(t) = sup H(t), 

x£C -Kr<! 

let us suppose that our statement is wrong, i.e., z = is a singular point. 
Then there are sequences Xk G C and — 1 < t k < 0, having the following 
properties: 

h(tk) = H(tk) = Mk = \v(xk,tk)\ —* +oo as k —> +oo. 

We scale our functions v and q so that scaled functions possess axial symme- 
try: 

u h (y, s) = X k v(X k y f , x 3k + X k y 3 , t k + \\s), X k = —, 
P k (y, s) = X\q(X k y\ x 3k + X k y 3 , t k + X 2 k s). 
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These functions satisfy the Navier-Stokes equations in Q(M^). Moreover, 

|uV*>0,0)| = l, y' k = M k x' k . (4.3) 

According to (j4.2p . 



i\ < M 
for all fceN. Thus, without loss of generality, we may assume that 

y'k^y* as k -► +°°- ( 4 - 4 ) 

Now, let us see what happens as fc — > +00. By the identity 

sup \u k (e)\ = 1 (4.5) 

e=(y,s)eC(M k ) 

and by (j4.ip . we can select subsequences (still denote as the entire sequence) 
such that 

u k ^u in L 00 (Q(a)), (4.6) 

and 

p k ^P in L|(Q(o)) (4.7) 

for any a > 0. Functions it and p are defined on Q- = K 3 x] — oo,0[. 
Obviously, they possess the following properties: 

ess sup \u(e)\ < 1, (4.8) 

sup (A(0, r; w) + £(0, r; u) + C(0, r; u) + D(0, r; p)) < A 1; (4.9) 

0<r<+oo >■ ' 

ess sup \y'\\u(y,s)\ < A 2 . (4-10) 

e=(y,s)eQ- 

Now, our aim is to show that u and p satisfy the Navier-Stokes equations 
Q_ and u is smooth enough to obey the identity 

1^,0,0)1 = 1. (4.11) 

To this end, we fix an arbitrary positive number a > and consider numbers 
k so big that a < M k /4. We know that u k satisfies the nonhomogeneous heat 
equation of the form 

d t u k - Au k = -divF k in Q(4a), 
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where F k = u k <8> u k + p k I and 

||^ fe |||,Q(4a) <ci(o) < oo. 
This is implies the following fact, see [13] . 

||VM fc ||| iQ(3a) <c 2 (a) < oo. 

Now, we can interpret the pair u k and p k as a solution to the nonhomogeneous 
Stokes system 

d t u k -Au k + Vp k = f k , divu k = in Q(3o), (4.12) 

where / fc = — w fc • Vu fc is the right hand side having the property 

ll/1f,Q(3a) <C 2 (o). 

Then, according to the local regularity theory for the Stokes system, see 
we can state that 

ll<VHf, Q{2a) + || W|| fiQ(2a) + ||VA; fc ||| iQ(2a) < c 3 (a). 

The latter, together with the embedding theorem, implies 

ll VMfc H3,f,(Q(2a)) + lb fc H3,f,Q(2a) < c ±( a ) ■ 

In turn, this improves integrability of the right hand side in (j4.12p 

ll/ fc H3,|,Q(2a) < c 4(0)- 

Therefore, by the local regularity theory, 

ll<VH 3 ,§,Q(2a) + II W|| 3) | iQ(2a) + ||VA; fe || 3i | iQ(2) < c 5 (a). 
Applying the imbedding theorem once more, we find 

l|V?/'|| 6i 3,Q(2a) + lb fc || 6 ,|,Q(2a) < ^{a) . 

The local regularity theory leads then to the estimate 

H^*ll6,|,Q(a) + II W|| 6i | iQ(a) + ||Vp fc || 6ifiQ(a) < C 7 (a). 
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By the embedding theorem, sequence u k is uniformly bounded in the parabo- 
lic Holder space C*(Q(a/2)). Hence, without loss of generality, one may 
assume that 

u k ^u in C*(Q(a/2)). 

This means that the pair u and p obeys the Navier- Stokes system and (14.111) 
holds. So, the function u is the so-called bounded ancient solution to the 
Navier-Stokes system which is, in addition, axially symmetric and satisfies 
the decay estimate ( 14.101) . As it was shown in [10] . such a solution must be 
identically zero. But this contradicts (j4.1ip . Theorems (13. ip and (13. 2p are 
proved. 



5 Appendix I: Proof of Theorem 12.4 

In what follows, we need a few known regularity results. 

Lemma 5.1. Assume that functions f G L m (B(2)) andq G L m (B(2)) satisfy 
the equation 

Aq = -divf m B{2). 

Then 

J \Vq\ m dx<c(m,n)( J \f\ m dx+ J \q - [q] B{2) \ m dx) . 

B{1) B(2) B(2) 

Lemma 5.2. Assume that functions f G L m {Q{2)) and u G W^°(Q(2)) 
satisfy the equation 

d t u — Au = f in Q(2). 

Then u G W^iQil)) and the following estimate is valid: 

\\dtU\\ m ,Q(l) + ||V 2 u|| miQ(1) < C(m,n) ||/|| m ,Q(2) + IMI^°(Q(2)) 

Proof of Lemmata 15. II and 15.21 can be found, for example, in p3j and ,_ 
Proof of Theorem 12.41 Step 1. Energy estimate. Take an arbi- 
trary number to < and fix it. Let k £ (z) be a standard smoothing kernel 
and let 

F £ (z)= J k £ (z-z')F(z')dz\ F = u®u, 
Q- 
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u ^z) = / k £ (z - z')u(z')dz' . 
Q- 

Assume that w E C^iQ-), where Q*° = R"x] - oo,£ [- Obviously, 

o 

w e G Co° (Q-) for sufficiently small e. Using known properties of smoothing 
kernel and Definition 12. 3[ we find 

J w ■ (d t u £ + div F £ - Au £ )dz = 0, VwG CffiQ*). 
Q- 

There exists a smooth function p e with the following property 

d t u £ + div F £ - Au £ = -Vp e , div u £ = (5.1) 

in Q°. Splitting pressure p e into two parts 

Pe=PF-+Pe- (5.2) 

and observing that the function Vpi?e is bounded in Q^, one can conclude 
that, by (IO) and (152)1 . 

Ap £ = in Q*_°, Vp e GL oo (g i0 ;R™). 

According to Liouville's theorem for harmonic functions, 

Vp £ (x, t) = a e (t), x e R n , -oo < t < to- 

So, we have 

<9 4 w £ + div F £ - Aw £ = -Vpjw - a e , div u £ = (5.3) 

in Ql°. 

Now, let us introduce new auxiliary functions 

t 
bet (t) = / a £ (r)rfr, t - 1 < t < t Q , 

to-l 

v £ (x,t) =u £ (x,t)+b £t0 (t), z= (x,t) EQ to . 



29 



Using them, one may reduce system (15.31) to the form 

d t v £ - Av £ = -div F £ - Vp F e , div v £ = (5.4) 

inQ*°. 

Let ip xo (x) = <p( x ~ x o) f° r a fixed cut-off function (p satisfying the con- 
ditions 

< ip < 1, tp = l in B(l), suppy? C B(2). 

To derive the energy identity, let us multiply (15. 4p by (f 2 v £ and integrate 
the product by parts. As a result, we have 

t 

2 |V7„. |2 t j./ 



A*) =: / ^ ( x )l^( a; ' t )r rf:r + 2 / I ¥ X0 \Vv £ \ 2 dxdt 

to-i 




/ 



f <fl o (x)\v E (x,t - l)\ 2 dx + f f A(pl o \v E \ 2 dxdt'+ 



t -l M n 



+ 



/ / (Pf« - tpjw]B(«o,2))u e • Vip 2 Xo dxdt'+ 



t -l 

t 



+ 

to- J 
Introducing 



/ /(^ £ " [F%(*o,2)) ■ -V&l^dxdt'. 



a £ (t) = sup / \v e (x, t)\ 2 dx 
x m n J 



B(x ,l) 

and taking into account that v £ (-,t — 1) = u e (-,t — 1) and |ii e (-, t — 1)| < 1, 
we can estimate the right hand side of the energy identity in the following 

way 

t 

I(t) < c{n) + c(n) / a £ {t')dt'+ 
to-i 

to t 

+ c ( ri )( / / l^ s - [pH B(a:o,2) 1 2 ^dtJ 2 f / a e (t')^j 2 + 

to-lB(x ,2) to-1 
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to . t 



+ c(n)( / / \F e -[F e ] BM \ 2 dxdty( [ f ip 2 X0 \Vv e \ 2 dxdt'+ (5.5) 



to-lB(x 0) 2) to-lR™ 

t 

+ / a £ (t')dt'Y, t - 1 < t < t . 



to-l 



Next, since |-F e | < c(n), we find two estimates 



to 

f f \F £ -[F% (x(h2) \ 2 dxdt<c(n) 



to-lB(x ,2) 

and 

to 

D (-oo,t ;-BMO(R n )) 



to 

/ / Ip** - [PF^B(x ,2)\ 2 dxdt < c(n)\\p F e\\ 2 Lc 



to-lB(x ,2) 



<cH||Fl^ (Q!o) < c (n). 



The latter estimates, together with (15.51) . imply the inequalities 

t 
a £ {t) < c{n) (l+ J a £ {t')dt^ , t - 1 < t < t 

*o-l 
and 

to to 

sup / / \Vv £ \ 2 dxdt < c(n) ( 1 + / a e (t)<it). 

z €R n J J ^ J ' 

to~lB(x ,l) t -l 

Applying known arguments, we can conclude 

to 

sup a £ (t) + sup / / \Vu e \ 2 dxdt < c{n). (5.6) 

t () -i<t<t x eR n J J 

to-l B(x ,l) 

It should be emphasized that the right hand size in (15.61) is independent of 
to- In particular, estimate (I5.6P allows to show 

sup b £to (t) < c(n). 

t -l<t<t 

31 



Now, let us see what happens if e — ^ 0. Selecting a subsequence if neces- 
sary and taking the limit as e — > 0, we state that: 

the estimate 

to 

ll&tollw*o-i,*o) + SU P / / \Vu\ 2 dxdt < c(n) < +00 (5.7) 
x eR n J J 

t -lB(x ,l) 

is valid for all t < 0; the system 

d t u to + divu ® u — Au = —Vp u ® u , divw = 

holds in Q to in the sense of distributions. 

The case to — can be treated by passing to the limit as to - ► 0. 
Step 2, Bootstrap Arguments By (15.71) . 

/ = divF = u ■ Vu G C 2 (Q-; M n ). 



Then Lemma 15. 1[ together with shifts, shows that 

Vp^e£ 2 (Q_;M n ). 
Next, obviously, the function u to satisfies the system of equations 

<V° - Au t0 = -u-Vu- Vp u ® u G £ 2 (Q_; W 1 ), 
which allows us to apply Lemma 15.21 and conclude that 

u t0 G W^{Q(zo, r 2 ); M n ), 1/2 < r 2 < n = 1. 
Moreover, the estimate 

iKl^Q^))^^,^) 

holds for any z = (x , to) 5 where x G M n and t < 0. Applying the parabolic 
embedding theorem, see [13], we can state that 

V^ = V«GW^(Q(* ,r 2 );R n ), 
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where 

1 1 1 

m 2 nil t^ + 2' 

By Lemma [5TTI by shifts, and by scaling, 



mi — 2. 



\V Pmu (;t)rdx<c(n,T 2 ,T^[ I \Vll(;t)\ m2 dx+l 
B(x ,r^) B(x ,r^) 

for 1/2 < Tg < r 2 . In turn, Lemma [5.21 provides two statements: 

u^ E W£(Q(zo, r 3 ); M n ), 1/2 < r 3 < r^ 

and 



|w*°n,^ 



»% a 1 «(*.-n.))^ c ( n ' T 3' 7 3)- 
Then, again, by the embedding theorem, we find 

Vu t0 = VueW^(Q(z ,T 3 );R n ) 

with 



1 



1 



1 



777,3 m 2 n + 2 

Now, let us take an arbitrary large number m > 2 and fix it. Find a as 
an unique solution to the equation 



a 



m 2 n + 2 

Next, for fc = [a] + 1, where [a] is the entire part of the number a, determine 
the number m fco+1 satisfying the identity 



1 



1 



An 



777 fco+ i 2 77 + 2 

Obviously, mk +i > m. Setting 

1 1 



T fc+ i = T fc - 



4 2 fc ' 



Tl = l, 



k — 1,2,,, 



and repeating our previous arguments k times, we conclude that: 

u t0 eW% o+1 (Q(z ,T ko+ i);R n ) 
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and 

l|w to || w 2,i lnl ,, < c(n,m). 

Since r^ > 1/2 for any natural numbers k, we complete the proof of Theorem 
12.41 Theorem 12.41 is proved. 

We can exclude the pressure field completely by considering the equations 
for vorticity uj = V A u. In dimensions three, differentiability properties of uj 
are described by the following theorem. 

Lemma 5.3. Let u be an arbitrary bounded ancient solution. For any m > 1, 
we have 

w = VAu6^ 1 (Q_;R 3 ) 

and 

dfUJ + u ■ Vcu — Auj = uj ■ Vu a.e. in Q_. 

Remark 5.4. We could continue investigations of regularity for solutions 
to the vorticity equations further and it would be a good exercise. However, 
regularity results stated in Theorem \5.3\ are sufficient for our purposes. 

Remark 5.5. Functions uj and Vuj are Holder continuous in Q_ and their 
norms in Holder spaces are uniformly bounded there, see [13]. 

Proof of Lemma 15.31 Let us consider the case n = 3. The case n = 2 
is in fact easier. So, we have 

dtOJ — Aco> = uj ■ Vu — u ■ Vuj = /. 



Take an arbitrary number m > 2 and fix it. By Theorem 12. 4[ 

|/| < c(n)(\V 2 u\ + \Vu\ 2 ) G L m (Q(z ,2)) 

and the norm of / in L m (Q(zo, 2)) is dominated by a constant depending only 
on m. It remains to apply Lemma T5.2I and complete the proof of Lemma 15731 
Lemma 15.31 is proved. 



6 Appendix II: Proof of Lemma 13.3 



According to the local regularity theory of the Navier-Stokes equations, see, 
for instance, [23], one can easily show that the pair v and q has the following 
differentiability properties: 

veW 2 /(C(a)x}-a 2 J -b 2 [), q G Wl'°(C(a)x] - a 2 , -b 2 [) (6.1) 

2 2 
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and 

v is Holder continuous in the completion of C(a)x] — a 2 , —b 2 [ (6.2) 

for any < b < a < 1. 

Now, we fix a number m > 2, multiply the equation for the velocity 
component u^ by ru\u\ m ~ 2 , where u = rv^, and integrate the product by 
parts. In view of (16.11) and (16. 2p . we find the following identity for u = \u\ m 

tfj 2 (x,t*)\uj(x,t*)\ 2 dx+ ^ m ~ - / 4> 2 \Vu}(x,t)\ 2 dxdt 
C -1 c 

u 



u{x, t)\ 2 (d t tp 2 + v ■ V^ 2 + A 2 ^ 2 + — ) dxdt. (6.3) 

-1 c 

It is valid for all —1 < £* < and for all cut-off functions vanishing in a 

neighborhood of the boundary of the space-time cylinder Cx] — 1, 1[. Here, 

A 2 ^ 2 = ip 2 rr+ip%3- So, (16.31) means that the energy norm of ipu is finite, i.e., 

u 

l^l^cxl-i.tJ — ess SU P / |V ;LL '( a; )^)| 2 ^ a; + / / \V{ipu)\ 2 dxdt 
te}-i,tJJ J J 

c -1 c 

u 2 

<c f f \u(x,t)\ 2 fd t ^ 2 + v -V^j 2 + A 2 ^j 2 + —^ + \V^\ 2 )dxdt (6.4) 
-i c 
for any — 1 < t* < 0. 

No, let us specify our cut-off function ip setting i/)(x,t) = $(x)x(t) and 
assuming that new smooth functions < $ < 1 and < x < 1 meet the 
following properties: 

supp$ E C{r\), $ = 1 in C(r), 

|V$| < — - — , |V 2( I>| < , \dtx\ < . 

r\ — r [r\ — r) z (ri — r) z 

Here, arbitrary fixed number r and r\ satisfy the condition 

13 , N 

-<r<n<-. (6.5) 
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If we let Q = C(ri)x] — rf, £*[, then 



\^\l£<-JL^I\u,\>dz + cI. (0.0) 



where 

/= / |^o;||a;|| : u|d2;. 

ri — r 



By Holder's inequality, 

J- / / , i 10 , \ 10 / / , , 5 , \ 5 / / i , i 10 \ 10 



/<— -( iM-rdz) f /H»cfe)-( /|^l3 



The left hand side of ( 16.51) can be evaluated from below with the help of the 
well-known multiplicative inequality 

\il)uj\~^dz\ <c\ipuj\ 2 Q. (6.7) 



Now, taking into account restriction (16. 5p . it is not so difficult to derive from 
( 16.61) and (16. 7p the following estimate 



3 rzr . /• .2 



\ijjuj\fdz) 10 <cM^-( / \Lu\ldzY. (6.8) 

Q 
Setting 



ri — r 

Q 



m = mk= {£) ' r i = rW = 2 + 2^+1' r = r(fc+1) 5 

V> = A, Qk = C(r W )x] - (r«)V*[, keN. 

one can reduce (16. 8p to the form 



(/|, H ^)^ < cM^il_(/ H ^)^ (6.9) 
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The only difference with respect to the usual Moser's technique is that one 
should take the limit as i* — > step-by-step. For example, for k = 1, the 
integral 

\u\ 3 dz 



Q(3/4) 

is finite and therefore we can pass to the limit as t* —>■ in (16.81) . Then 
we may pass to the limit as t* — > in (16.81) for k = 2 and so on. Tending 
k —>■ +oo, we complete the proof of Lemma I3TB1 in more or less standard way. 
Lemma 13.31 is proved. 
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